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We demonstrated that the chiral soliton lattice formed out of a chiral helimagnet exhibits coherent 
sliding motion by applying a time-dependent magnetic field parallel to the helical axis, in addition 
to a static field perpendicular to the helical axis. To describe the coherent sliding, we use the 
collective coordinate method and numerical analysis. We also show that the time-dependent sliding 
velocity causes a time-varying Berry cap which causes the spin-motive-force. A salient feature of 
the chiral soliton lattice is appearance of the strongly amplified spin motive force which is directly 
proportional to the macroscopic number of solitons (magnetic kinks). 



I. INTRODUCTION 

Spin-based-electronics (spintronics) is now an emerg- 
ing field. An essential notion behind this emergence is 
the fact that the 'spin magnetic current' without rely- 
ing on electric current would magnificently reduce the 
energy loss and switching time during the information 
read/write processes. At the heart of spintronics is to 
drive motion of magnetic textures in a controllable man- 
ner. There are two ways to drive the motion, i.e., inco- 
herent and coherent methods.^ The incoherent method 
is typically realized by injecting a spin-polarized cur- 
rent into a sample.^ On the other hand, the coherent 
method is realized in a magnetically ordered state by 
twisting the phase angle of the magnetic order parame- 
ter which directly couples to a magnetic field. In the co- 
herent method, the phase rigidity (stiffness) of the whole 
spin system makes it possible to transmit the spin ro- 
tation at one end of a sample to the other end via spin 
torque transfer. Although the coherent method has po- 
tential advantage, it has not been extensively studied in 
the context of the present-day spintronics, because it is 
not so easy to prepare rigid phase object which transports 
experimentally measured quantity in a controllable man- 
ner. To seek for the possibility of the coherent method is 
interesting. 

In this paper we propose that chiral helimagnets are 
promising candidates to realize the coherent method. 
The chiral helimagnetic (CHM) state is character- 
ized by the vector spin chirality as an order parame- 
ter. The structure is stabilized by the antisymmetric 
Dzyaloshinskii-Moriya (DM) interaction and realized in 
crystals without rotoinversion symmetry. A guiding prin- 
ciple for materializing this effect is symmetry-adapted 
material synthesis, i.e., the interplay of crystallographic 
and magnetic chirality plays a key role there. 

The CHM state is, however, nothing more than non- 
colinear linear (harmonic) spin structure. To carry phys- 
ically measurable quantity, we need nonlinear structure. 



Fortunately, under the magnetic field applied perpendic- 
ular to the helical axis, the CHM is transformed to a 
non-linear magnetic structure called a chiral soliton lat- 
tice (CSL) [see Fig. [1] which is equivalent to a magnetic 
kink crystal (MKC).'^'^ In the CSL state, the ground state 
possesses a periodic array of the commensurate (C) and 
incommensurate (IC) domains partitioned by discom- 
mensurations. Recently, using Lorenz microscopy and 
small-angle electron diffraction, the CSL was experimen- 
tally verified in a hexagonal chiral magnet Cri/3NbS2i^ 
Present authors have discussed physical outcome of the 
CSL state from various viewpoints.—"— 

As pointed out in Refs.— i^, the CSL exhibits coherent 
collective sliding motion in non-equilibrium state. Once 
the sliding is triggered, the CSL maintains its persistent 
motion assisted by the dynamical generation of the iner- 
tial mass. The mass generation is understood by Doring- 
Becker-Kittel mechanism of the moving domain wallr^"— 
In this mechanism, the longitudinal (out-of-plane) com- 
ponent of the slanted magnetic moment inside the do- 
main wall emerges as a consequence of translational mo- 
tion. An additional magnetic energy associated with the 
resultant demagnetization field is interpreted as the ki- 
netic energy of the wall. 

The incoherent current injection method to drive the 




FIG. 1: Schematic picture of the CSL. SMF generation needs 
static transverse field Hx and time-dependent longitudinal 
field Hx. 
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sliding was already proposed by the present authors.^ In 
this paper, we demonstrate that crossed magnetic fields 
are eligible to cause the coherent motion of the whole 
CSL. Here, we mean by the crossed fields that in addition 
to a static field perpendicular to the helical axis, which 
stabilizes the CSL formation, a magnetic field parallel to 
the helical axis is imposed. 

Once the coherent motion occurs, the natural question 
arises as to whether the motion has observable conse- 
quences for the spin motive force (SMF).— Quite natu- 
rally to expect that the time dependence of the longitu- 
dinal magnetic field manifests itself in a temporal regime 
of the SMF. Time dependences of the spinmotive forces 
are classified into three types, i.e., (i) transient, (ii) con- 
tinous ac, (iii) and continous dc ones. For example, the 
domain wall motion^^ and electron transport through fer- 
romagnetic nanoparticlesi^ lead to the type (i) SMF. A 
vortex core dynamics of a magnetic disk caused by an 
oscillating magnetic field directed in the disk plane in- 
duces a continous ac spinmotive force of type (ii)-'^^ A 
resonant microwave excitation of a comb-shaped ferro- 
magnetic thin film produces a continous dc spinmotive 
force of type (iii)'^ We will demonstrate that the time- 
dependent longitudinal field, as shown in Fig.[TJ possibly 
causes the SMF of type (i) and (ii) in the chiral helimag- 
net. As a remarkable feature, we note that our CSL is 
a macroscopically ordered object, which contains macro- 
scopic amounts of magnetic solitons. Due to this huge 
number of solitons, the SMF is expected to be strongly 
amplified as compared with the SMF caused by a single 
magnetic domain wall in a ferromagnet. 

In Sec. I, we present a model and summarize necessary 
background on the CSL dynamics. In Sec. II, we demon- 
strate that under the presence of the longitudinal field, 
the CSL becomes unstable and coherent motion occurs. 
In Sec. Ill, we present numerical analysis of dynamics to 
support analytical consideration presented in Sec. III. In 
Sec. IV, we discuss the SMF associated with the coherent 
motion. We conclude our results in Sec. V. 



II. STATIC DEFORMATION OF CSL UNDER 
CROSSED MAGNETIC FIELDS 

A. Basic equations of chiral soliton lattice 

1. Static structure 

Mono-axial chiral helimagnet is described by an effec- 
tive one-dimensional Hamiltonian, 



H 



-JJ2S^■S,+ l-D■J2 S^ X S^+,+H-J2 (1) 



where Si is the local spin moment at the site i, J > is 
the nearest-neighbor ferromagnetic exchange interaction, 
D = Dez is the mono-axial Dzyaloshinskii-Moriya (DM) 
interaction along a certain crystallographic chiral axis 



(taken as the 2;-axis). We take z-axis as the mono-axis 
and apply magnetic field H — gfisH — g^B{Hx, 0, Hz) 
in the xz-plane, where g is the electron g-factor and 
Ms — ti/2m is the Bohr magneton. 

In the semi-classical approach, because of the slowly 
varying nature of the spin variables, it is legitimate to in- 
troduce the continuous field variable S{z) = Si5{z — 
Zi) = Sn{z). A unit vector field 'n{z) is represented as 

n(z) — [sm9{z) cos (p{z), sin 9 (z) sin (p{z), cos 9 (z)], (2) 

by using the polar angles 9{z) and (p{z). The continuum 
version of the Hamiltonian H — J^^ dzH , where L 
denotes the linear dimension of the system, includes the 
Hamiltonian density. 



n 



Zan 



D ■ n X dzTb 



^H-n. (3) 



Here ao is the atomic lattice constant along the chiral 
axis (ao — lOA in Cri/3NbS2.^) 

For H = 0, the Hamiltonian ([T|) gives a a;j/-planer he- 
limagnetic structure, 9{z) ~ 7r/2 and (p{z) = Q^z with 
the modulation wave number being given by Qq ~ 
ag ^ arctan(£'/J) ~ a^^D/J. 

For only non-zero transverse field perpendicular to the 
chiral axis, the CSL structure becomes the ground state 
characterized by = 7r/2 and 

fo{z) — 2am ( z 



where am is the Jacobi's amplitude function with the 
elliptic modulus k (0 < k < 1), and E = E{k) is the 
complete elliptic integral of the second kind. The elliptic 
modulus K is determined from the condition, 




AE 
nQoao 



This equation is also written as Hx/Hc ~ [i^/E) 
introducing the critical field corresponding to k = 1, 



Hr 



(4) 

by 

(5) 



at which a incommensurate-to-commensurate phase 
transition occurs. The spatial period of the CSL is given 
by icsL = S/fE'/TrQQ, which continuously increases from 
^CHM =27r/QQ to infinity when the magnetic field in- 
creases from zero to H^- Here, Lchm is the spatial pe- 
riod of CHM under zero field. K = K{n) is the complete 
elliptic integral of the first kind. 



2. Individual spin dynamics 

Next, we write down the basic equations for dynam- 
ics. Using the Hamiltonian density ([3]), we construct the 
Lagrangian density. 



C = (cos9 — 1) dtf ~ H. 



(6) 
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To incorporate the damping effect, we use the Rayleigh 
dissipation described by 



and 



(7) 



with a being a smah coefficient specifying the Gilbert 
damping. The Euler-Lagrange equations of motion is 
then given by 



hS . SH SW 
smOdtO = h -— , 

Qq dip 0(p 

hS . SH SW 



(8a) 
(8b) 



which leads to the Landau-Lifshitz-Gilbert (LLG) equa- 
tions, 

hSsinddtO = -JS^al {sin^ 9dlip+ s.in2e {d,ip) {8,9)] 
+ DS'^ao sin 20 (8,6) -H^Ssin 9 sin ip 
+ ahS sin^ e (dtip) , (9a) 

+ DS^act sin 29 {8ztp) — HxScos 9 cos (p 

+ H,Ssm9-ahS8t9. (9b) 

The LGG equations describe the individual (not collec- 
tive) spin dynamics. 

3. Collective dynamics 

To consider the sliding motion of the CSL, we use the 
collective coordinate methodic used introduced in Ref4. 
In this formulation, the CSL dynamics is fully described 
by two dynamical variables, the center of mass position, 
Z, and the out-of-plane quasi-zero mode (OPQZ) coor- 
dinate ^0- Using them the sliding solution is written as 



ip{z,t) = ipo[z- Z{t)] , 
9{z,t)^TT/2 + ^oit)uo [z-Z{t)] 
The zero-mode wave function 



(10a) 
(10b) 



uo{z) 




T^Qo 



KE 



(11) 



serves as the basis function of the ^-fluctuations localized 
around each soliton and ^o(0 the OPQZ coordinate. 
Here dn is the Jacobi-dn function. The function uq(z) 
exactly corresponds to the topological charge distribution 
because 8zPo{z). Using these variables, the Lagrangian, 
£ = J dzCdz, and the Rayleigh term, W = / dzW, are 
respectively written as 



hS ■ e'i^ 
C = — 3-/C^o^ rCo 



(12) 



ahS 
'24 



(13) 



where JC — dzuo{z)8z'Po{z) and K, = dz {8z(po{z))'^ 

as given m Ref4. Furthermore, e^^^ ~ D'^S^/2J IS an 
energy gap of the 0-mode caused by the DM interaction 
which plays a role of easy-plane anisotropy. We here note 
a useful relation 



(14) 



[see Eq. (16) of Ref.^] 

Under the condition of weak field, k <C 1, we have 
K. ~ QoVl and M ~ Q^L. We have ^o{t) + only for 
nonequilbrium state where the CSL exhibits sliding mo- 
tion. An emergence of such coherent collective transport 
in non-equilibrium state is a manifestation of the dynam- 
ical off-diagonal long range order. Using Eqs. (jlOap and 
(jlObl) . Eqs. ([^ and (|9b)) lead to equations of motion for 
collective dynamics, 

?i/C^o = -ahMZ, (15a) 
hK.Z = 2S-^ei^'^^o + ahia - VlHz. (15b) 

This set of EOMs differs from Eqs. (37) given in Ref.* in 
that the sd interaction is absent and the longitudinal field 
is present. In Ref.®, the incoherent driving of the CSL 
sliding motion was driven by the spin-torque-transfer 
from the spin-polarized current to the local spins. On 
the other hand, in the present case, we are discussing the 
coherent driving caused by the uniform time-dependent 
magnetic field Hz- 

Eqs. (|15al) and (I15bp are readily solved to give 



Z{t) = Ce-*/^c^L 

p — t/ ''"CSL 



^ (l + a2)/i 



^t'/rcs.^E^dt', (16a) 
dt' ' V y 



e^'/''cs^Hz(t')dt', 



(16b) 



where we used the relation M^jK? — 1 for a weak field. 
Constants C and D are determined by the initial condi- 
tion Z(0) = and ^o(O) = 0. The intrinsic relaxation 
time of the CSL, caused by the Gilbert damping, is in- 
troduced by 



TCSL 



2e 



S 



which is also written as 



I/tcsl 



(17) 



(18) 



with uja 



gap — Eq /h being a characteristic frequency of 
the gap. In the case of static Hz, we have trivial relax- 
ational dynamics where the sliding motion never persists. 



Furthermore, the DM interaction D gives rise to a finite 
relaxation time. Eq. (|16bp means that the longitudinal 
field first directly coupes to and drives its grows via 
the Gilbert damping process. Then, the sliding motion 
follows the growth of ^o- This process is consistent with 
intuitive ideas developed by D6ring.<i2' 

We here emphasize that the two coordinates and 
Z coupled to each other via the Gilbert damping a [see 
Eq. (|15a|) ]. If there were no damping, we would have no 
correlated dynamics. This fact means that the CSL never 
realizes dissipationless collective motion. As we will see 
in Sec. V and appendix C, the presence of the damping 
is essential to drive SMF. 



e{z) 

l.b/uf 
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FIG. 2: Spatial modulation oie{z) = 7r/2+6'i(2) for H^/H^ = 
0, 0.1, 0.2, 0.3 (from above to below) with keeping H^/ = 
0.1. Here a dimensionless coordinate z = nQoz/AE is used. 



4- Comparison of material parameters to theoretical 
formulae 

It may be worthwhile to summarize relation between 
experimental data and theoretical parameters by taking 
an example of Cri/3NbS2^ and theoretical formulae. In 
this sample, it is reported that ap — 1.212 x lO^^m, 
icHM = 27r/(5o = 4.8 X 10~^m. Cr atoms are in the 
trivalent state and have localized electrons with spins of 
5 = 3/2. An observed critical field is He = 2300Oe corre- 
sponding to 0.31K. We have an estimation for the ratio, 
D/J — tan((5oao) = 0.16. Another important quantity 
is the excitation gap which is estimated as eQ^-* ~0.38K by 
using Eq. ([T4| . The intrinsic relaxation time of the CSL 
is also estimated as rcsL — {a + a~^) x 3.0 x 10~^^sec. 
A small damping such as a ~ 10~^ leads to tqsl ~ 
3.0 X lO^^sec. Smaller damping causes longer relaxation 
time. To realize a longer period of relaxation processes, 
it is desirable to realize a smaller value of a and a smaller 
gap frequency i:jgap. 



B. Static deformation of CSL 

Now that we have set up all the necessary equations 
for dynamics, we proceed with the stability analysis of 
the CSL against crossed magnetic fields. Before going 
to dynamical deformation, it is worth while to study the 
static deformation. For analysis, we here consider the 



weak field limit. 



< 1. 



< JS and 

We introduce the static deformations as 6{z) = tt/2 + 
s6(z) and (/^(z) = 1^90(2^) + s(p{z), with the small fluctu- 
ations 0, (p and s being a dummy expansion parameter. 
As derived in appendix A, we have 



(19) 



where the Wronskian W is given in Eq. (jA7p and the 
coefficients Ci^2 are given by Eqs. (lAlip and (jA12|) . The 
functions <^i.2 are a pair of the fundamental solutions of 



Lame equation [Eq. (|A2[) ] . In Fig. [5] we show the ob- 
tained spatial modulation of 6{z) = 7r/2 + 9i{z) for vari- 
ous Hz /He with keeping Hx/Hz = 0.1. Since the coeffi- 
cients ci.2n fall exponentially with growth of n, we retain 
only terms with n = ±1 which dominate the terms with 
|n| > 2. We see that finite Hz tends to orient the spins 
toward the z-direction but causes non-uniform spatial 
oscillation. This oscillation implies the static defomation 
considered here is unstable against dynamic deformation. 
We will see this dynamic deformation corresponds to the 
sfiding motion of the CSL. 



III. COHERENT SLIDING DYNAMICS BY 
CROSSED MAGNETIC FIELDS 

A. Energy and momeutum associated with the 
sliding motion 

In this section we will show that the static deformation 
of the CSL is unstable against the dynamical instability, 
i.e., spontaneous coherent sliding motion of the whole 
CSL. As the present authors previously pointed out^^ 
the Lagrangian constructed from the Hamiltonian Eq. 
^ has the hidden Galilean symmetry induced thorugh 
Lie analysis. 

This hidden symmetry justfies that the CSL has a lin- 
ear momentum, 



Pz ^ — I dz{l- cos 6*) dzip, 



(20) 



associated with the kinematic Berry phase. We connect 
the momentum variation SPz to the energy variation)^! 



SE 



(21) 



associated with the sliding motion. In the absence of the 
dissipation, plugging the EOMs (|5al) and ((8b)) into Eq. 
m, we have 



6E^hS 



dz siii6 [ —Sip 

^ dt ^ dt 



(22) 
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The sliding motion means 9 and Lp are funcitons of m = 
z — Z[t) and we can make the replacement, dt — >■ —Zd^ 
dz ^ du- Then, we easily obtain the relation, 



6E = ZSP, 



(23) 



Based on this relation, we see that the coherent sliding 
occurs if the condition ZSPz < is sataified for a given 
momentum transfer 6Pz from the environment. We will 
see that in the present case the longitudinal magnetic 
field Hz gives the momentum transfer and drives the slid- 
ing motion. 



B. Coherent sliding caused by a transient 
longitudinal field 

We first consider a transient longitudinal field 



Hz{t)^Hzoil 



(24) 



switched on in addition to a perpendicular static field 
Hx which stabilizes the CSL. Eq. (I16ap gives the sliding 
velocity. 



Z{t) = -Vo ^^^^ fe-*/rcsL _ g-*/T 
TcsL -TV 

where the characteristic velocity is defined by 



Hzo 



hQoil 



(25) 



(26) 



Eq. (|25l) indicates Z < for T < tcsl and then the 
condition ZSPz < is satisfied. It is to be noted that 
if the chirality of a crystal were inverted, i.e., Qo is in- 
verted to — Qo: the velocity would be inverted. So, the 
sliding orientation and the crystal chirality correlate to 
each other. 

In the case of Cri/3NbS2, The characteristic velocity 
is estimated as Vq ~0.13[m-s~^-Oe~^]. So, the sudden 
switching of the longitudinal magnetic field Hz , satisfing 
the condition T < tcsl, will easily cause the coherent 
sliding motion the CSL. 

In Fig. [21 we show time-evolution of the sliding veloc- 
ity. We sec the velocity grows linearly with time shortly 
after the field Hz is switching on. Then, after the relax- 
ation time of the filed, T, the velocity begins to relax. 
Then, it finally relaxes to zero over the time scale of 
the Gilbert damping tcsl- Therefore, to realize longer- 
lasting sliding motion, a smaller value of a and a smaller 
gap frequency ujgap may be desirable. 



C. Coherent oscillating motion under AC field 

It is also possible that an oscillating longitudinal field 
Hz{t)^Hzisin{nt), (27) 




FIG. 3: Time dependence of (a) longitudinal field Hz{t) = 



H,o(l 



and (b) velocity Z /Va for T = 0.5 tcsl- 



causes a coherent oscillating motion of the CSL in addi- 
tion to a perpendicular static field Hx- In this case, Eq. 
Eq. (|16ap gives the velocity, 

Z{t) = Vi [e~*/"'°"^ - f^rcsL sin(m) - cos(f7i)] , (28) 
where the characteristic velocity is defined by 



Vi = 



Hzi^Tcsi. 



(29) 



Unlike the case of the transient filed, the oscillational slid- 
ing motion is sustained as a long-term stationary state. 
This is because in the AC case the energy associated 
with the CSL motion is perpetually supplied by the AC 
field. It is also seen that the Gilbert damping causes out- 
of-phase oscillation[cos(i7t)]. For experiment, it may be 
useful to note 



^0 



(1 



^T-QSL Hzl 



(30) 



'CSL 



In Fig. m we show oscillating response of the sliding 
velocity to the longitudinal AC field. The transient state 
rapidly relaxes over the time scale of T, to the stationary 
forced oscillation with a phase shift due to the damping. 

Here we comment on the relation of the sliding dynam- 
ics to the electron spin resonance (ESR).^ Provided the 
whole CSL is in a state of rest, the longitudinal AC field 
is able to excite the small amplitude phonon-like mode 
(chiral soliton lattice phonon or magnetic kink crystal 
phonon) propagating over the CSL. However, because the 
AC field is uniform, the resonant phonon absorption oc- 
curs only when the momentum absorbed by the phonon 
coincides with the reciprocal vector of the super-Brillouin 



6 



(a) Hzit) 



1 \2 




FIG. 4: Time dependence of (a) longitudinal field Hz{t) — 
Hzi sm{^lt) and (b) velocity Z /Vo for Q^^ = O.Stcsl- 



zone of the CSL. Once the resonance condition is satis- 
fied, the microwave energy would be consumed to excite 
the CSL phonons. On the other hand, in the case of off- 
resonance, the sliding motion would be driven. We also 
note that the excitations associated with the fluctuations 
of if is totally irrelevant to the CSL phonon excitation. 
On the other hand, the sliding motion is a consequence 
of the correlated dynamics of coupled 9 and (p. 



IV. NUMERICAL ANALYSIS OF DYNAMICS 



ft 



46000 46050 46100 46150 

FIG. 5: Numerical dependencies of 9 on the coordinate z. 
The length of the chain is L = 10^ sites. The calculation 
is performed with the same parameters as in Fig. (2] The 
perpendicular field, Hx/Hc — 0; 0.1; 0.2; 0.3 (from above to 
below), is normalized to the critical field He- 

where S = arctan(£)/ J), ~ H^/JS, and Pz = Hz/JS. 
Here, the time scale Tq = h/JS, and the dimensionless 
time T = t/TQ is introduced. 

First we consider static spin configurations. In order to 
perform numerical computations, we adjust the problem 
to the form convenient for an iterative routine, i.e., 

X + B,-if + (C,_i -I- P.+i + P^f]'^'^, 

(32) 

cos ipi = {Ci-i + Vi+i + 13^) 

X [(A+i + B,-, f + (C,_i + 2?,+i + /3.)']-'/', 

(33) 





A. Static deformation of CSL 

So far we discussed the CSL dynamics in analytic man- 
ner. To justify the obtained results, it is desirable to 
perform numerical simulations of the dynamics. For nu- 
merical analysis, we start with the lattice version of Eqs. 
(l9al) and (l9bl) written as 



= Y 1 + -j2- sme'i-i sm(v?i - -|- 6) 



- a/1 + -j2- sm^i+i sin((^i+i ~ (p^ + S) 

- sin 

^ = -(cos^j+i +cos6'i_i) 
ar 



/3x cot 9i cos (fit + /3z 



(31a) 



1 H ^ cot 6i sin Oi^i cos{(pi — (pi-i + S) 



1 + — J- cot 6i sin^i+i cos((/3j+i - (fi + S) 



(31b) 



cos di = (cos 9i+i + cos 6i^i + Pz) 
X [(cos Oi+i + cos 9i^i + jizf 

+ (A+i + s,_i)' + (c,_i + p.+i + pxfr^'\ 

(34) 



where we defined 




v/r~KD7J2sh 



f sin {(pi + 5) 

1- (35) 
cos {(pi — d) ' ^ ' 

\ cos [pi + 5) 



The spin configuration is found by using the original 
spin variables and iteratively re-pointing each along the 
effective local field due to its neighbors. Scanning linearly 
through the chain, spin variable at each site is updated 
in sequence, being reset along the net field due partly to 
some unchanged neighbors and some that have already 
been re-pointed. This gives convergence faster than a 
synchronized global update. 

The most difficult computational problem in carrying 
out this program is to find the initial configuration that 
relaxes to a target spin configuration. It is meaningful to 
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46 000 46 050 46100 46150 



FIG. 6: Coordinate dependence of the tp (shown as mod 2-k) 
obtained numerically for the chain of length L = 10"" sites at 
the same parameters as in Fig. [2] and Hx/Hc ~ 0.1. From 
the almost strict linear behavior, one sees that ip acquires 
almost no change. 



1.56907 



1.56906 




1.56905 .. 

1 .56904 \ ^» 

16 17 18 19 20 21 



FIG. 7: A comparison between numerical data for the static 
case (dotted line) and analytical expression given by Eq. (|19|) 
(solid line). The fields are fi^ — b cos St and j3z — b sin St, 
6=10"*, St = -K/18. 



(Fig.|6]). The numerical data are imposed on the theoret- 
ical predictions [Eq. (|19p ] as shovifn in Fig. [T] Evidently, 
they reveal a good agreement with each other (Fig. [7]). 



B. Dynamics 

A search for dynamical solutions is carried out by using 
the eighth-order Dormand-Prince method implemented 
in Ref.^^. The embedded Runge-Kutta integrator with 
an adaptive step-size control ensures a relative tolerance 
10"^^. The integration spans a period of time from zero 
till 2 X IO'^tq. The length of chains used in the com- 
putations amounts to 10^ -I- 1 sites. A time evolution 
of magnetization was monitored by recording the time 
dependencies of the 9 and (p variables for the central site. 

In Fig. [51 we show numerical result for the velocity. It 
is seen that the velocity decreased linearly with increasing 
(^z field. In the calculations (^x was held constant. 

Figs, ini [TU] show the time dependence of 9 and (p, re- 
spectively, under the oscillating for different values 
of the damping parameter a. These calculations make 
evident a salient feature of the forced oscillations. In 
the initial stage of time evolution, the longitudinal mag- 
netic field excites intrinsic eigenmodes that are superim- 
posed on the field-driven oscillations. The eigenmodes 
fade away in the steady-state regime and they damped 
more rapidly the greater the parameter a. The period 
of the forced oscillations 27r/Sl exactly corresponds the 
period of the driving force. In appendix B, we present 
the detailed analysis of this forced oscillation in line with 
the numerical analysis. 



impose appropriate boundary conditions too. Obviously 
this is a rich problem with a wide choice of options. In 
our simulations we choose a starting configuration as the 
simple spiral, ipi = q{i — zq) and 9i = 7r/2, and take the 
free boundary conditions. The coordinate zq corresponds 
to a position with tpi — n in a middle of the chain of 
the length L. The value of g = 0.1 is taken throughout 
the numerical calculations to make more appreciable a 
spatial modulation of solutions. 

We revealed that a convergence of the iteration process 
is very slow. The iterations stop if the sum 



\ 



(fc-i) 



L 

E 

1=1 



(36) 

taken over the chain on the fc-th step is less than tolerance 
10~*. To reach the accuracy, around 205 x 10^ iterations 
are required. 

The numerical behavior of 9 shown in Fig. [5] repro- 
duces the theoretical finding (see Fig. (2]), i.e. an increas- 
ing of the longitudinal field enhances a modulation of a 
conical structure along the z axis. The calculation con- 
firms another assumption of the Sec. IV, namely, the 
longitudinal field causes no changes in the variable (f 



V. SPIN MOTIVE FORCE 

Now that we have obtained the CSL dynamics under 
the crossed magnetic field, we will go on to discuss pos- 
sible SMF generation. Because the sliding motion of the 




FIG. 8: (a) Linear time dependence of pit) for the central 
site of a chain of length L = 1000001 obtained by numerical 
simulations. The fields are = b cos St and ~ mb sin St, 
b — 10"'^, 5b — it/ 18. With a growth of (or m) a slope 
of the curves increases linearly according to analytical result. 
The lines correspond to m = 1.0 — 1.4 with the step 0.1 from 
above to below, respectively, (b) The dependence Z on the 
factor m extracted from the ip{t) data. 
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FIG. 9: The time-dependent variation of the 6 for different 
damping parameters: (a) q = 0, (b) a — 0.01, (c) a — 0.02, 
(d) a — 0.05. The fields are fix = bcos^;, and jSz — —b sin St, 
b = 10"^, St = Tv/18, the ratio /J^/H equals to 0.01. 
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5000 10000 15000 2000O 



5000 10000 15000 2000O 



FIG. 10: The time-dependent variation of the intrinsic mode 
of ifi for different damping parameters: (a) a = 0, (b) a = 
0.01, (c) a = 0.02, (d) a = 0.05. The fields and frequency are 
the same as in the previous Figure. 



Then, we obtain a general expression for the SMF is given 

by 

ea{t) = dzE^{z,t) = (^J cosedifi^ , (38) 

where the contour T is taken on the sphere presenting a 
space of the order parameter n. The vohage is related 
via the Stokes theorem with a change of area (Berry cap) 
S on the sphere enclosed by the contour,— 



eait) 



ha dS 



(39) 



This involves an analogue of Faraday's law for the emer- 
gent electromagnetic field, where a magnetic field of a 
Dirac monopole with a charge h/2 plays a role of the flux 
enclosed by the Berry cap S. 

In the present case of the CSL dynamics, by using the 
collective representation [Eqs. (|10ap and (jlObp ]. the SMF 
for majority spins (the case of minority spins has the 
opposite sign) is computed as 



Ui){z)dz, 



(40) 



This expression indicates that the SMF arises only for 
^0 7^ 0, i.e., the time-dependence of the Berry cap is 
essential to cause the SMF [see Fig 111]. This observation 
is consistent with the discussion given in Ref.— . 
Using Eq. (ITT|) . we obtain 



UQ{z)dz = — 




KE 



L 



■Q, 



where 



ipo{L) - v?o(0) 
27r 



(41) 



(42) 



is the topological charge representing the number of soli- 
tons over the whole length of a sample. Using the rela- 
tion, ^0 = —aQoVLZ , obtained from Eq. (llSap . Eq. 
(l40l) finally reduces to 



CSL accompanies the dynamical deformation of the spin 
texture, we naturally expect the SMF to occur in the con- 
figuration presented in Fig. [TJ Generally speaking, when 
conduction electrons adiabatically see a spatially modu- 
lated spin structure along the z axis, spinor wave function 
locally follows the background. Consequently, the spinor 
space turns out to be curved. The corresponding curva- 
ture is represented by the gauge (Berry) connections. In 
the adiabatic picture, the Berry curvature in the spinor 
space acts as an effective electric field,— >2i 



ha 
Te 



ha 
27 



sin 9 [dzOdtip - dz(pdtO) 



(37) 



e{t) ~ Q-7TaQoZ{t), 

e 



(43) 



where we used the relation V KE ~ 7r/2 in the case of a 
weak transverse field. As expected, the SMF is directly 
proportional to the macroscopic number of soliton, Q. 
It is worth whole to compare the obtained formula [Eq. 
(l43t ] with the one used for the SMF induced by domain 
wall motion. •^^ In the present case of CSL, the SMF to 
be strongly amplified by the prefactor Q. Furthermore, 
it should be stressed that the dissipative dynamics is es- 
sential to drive SMF. Actually, the SMF is proportional 
to the Gilbert damping parameter, a. In appendix C, 
we show that the dissipationless rigid motion of the CSL 
never produces SMF. 
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(a) _ (b) 




FIG. 11: Switching of the time-dependent longitudinal field 
Hz{t) causes the change in Berry cap 5 from (a) to (b) in 
a time-dependent manner. This time-dependence causes the 
SMF along the chiral axis. Each allow represents local spin 
configuration n{z). The Berry cap is associated with the area 
traced out by n{z). 



In the case of the time-dependent longitudinal field, 
H,{t) = iJ^o(l-e~*/^), plugging Eq. ^ into Eq. O, 
we immediately obtain 

s{t) ^ ~So^^ (e-*/-- - e-*/^) , (44) 

TCSL - J V / 

where 

eo = aQ-nQoVo. (45) 
e 

In Cri/gNbSa, using Qo ~ 1.3 x 10^[m-i] and 
Vq ~0.13[m-s^^-Oe~^] as estimated in Sec. III.B and as- 
suming a ~ 10"'^, we have an estimation, Eq ~ 0.36QHzo 
[nV] when H^o is measured in Oersted. In the case where 
the sample size is L ~lmm, the upper bound of Q along 
the helical axis amounts to L/Lqum — 10^. Therefore, 
we expect that eo amounts to l[mV] for Hzo — lO^Oe as 
an example. To experimentally sustain the SMF, it may 
be desirable to apply a sequence of the pulse fields. 

Here we comment on the physical reason why the SMF 
is proportional to the Gilbert damping factor a in Eq. 

An essential point is that emergence of the SMF is 
a direct consequence of time-varying Berry cap,^^ which 
needs finite ^o- Now, as is clearly seen from basic EOMs, 
(|15al) and (|15bl) . if the Gilbert damping were absent, 
^0 and Z are dynamically decoupled and consequently 
£.0 — £o{0) — for all the time, i.e., the sliding motion 
can never be sustained. This situation is totally differ- 
ent from the case of 180-degree Bloch wall. The CSL 
is regarded as an array of 360-degree walls and we need 
some mechanism which enables the magnetic moments to 
rotate around the chiral axis. Only one possible mecha- 
nism to realize this rotation is the Gilbert damping pro- 
cess. This is the reason why the SMF is proportional to 
a. More intuitively speaking, at the first stage the lon- 
gitudinal field Hz directly couples to and causes the 
out-of-plane canting of the magnetic moments [see the 
Lagrangian (IT^ ]. At the second stage, because of the 



Gilbert damping, the magnetic moments start damped 
precession to relax back into their original directions. 
This motion triggers the moments to rotate around the 
chiral axis and eventually leads to the collective sliding. 



VI. CONCLUDING REMARKS 

In this paper, we demonstrated that the chiral soliton 
lattice (CSL) exhibits coherent sliding motion simply ap- 
plying a time- dependent magnetic longitudinal field, in 
addition to a static transverse field. The driving force 
of the sliding is given by the Zeeman coupling of the 
collective coordinate with the longitudinal field. This 
mechanism is intuitively understood by DBK mechanism 
of the moving domain wall (DW)i^^— In the DBK mech- 
anism, once the domain wall begins to move, so-called 
demagnetization field is dynamically generated inside the 
wall. The demagnetization field supplies spin torque to 
keep the inertial motion. Actually, the DBK mechanism 
of a single DW was analyzed by exactly the same proce- 
dure presented in this paper. In the CSL dynamics, the 
longitudinal magnetic field kicks off the demagnetization 
and drives coherent sliding motion. To demonstrate the 
coherent sliding motion, we first used the collective co- 
ordinate method, and then confirmed the result by com- 
putational analysis. 

The time duration of the sliding motion is charac- 
terized by the intrinsic relaxation time determined by 
Eq. (1181) . Providing the gap energy e'^^ varies between 
0.1 ~ lOK and the Gilbert damping constant varies be- 
tween lO"'^ 10^2 

we expect the time duration varies 
between 10~^ ^ 10~^ sec. To realize longer-lasting slid- 
ing motion, a smaller value of a and a smaller gap fre- 
quency cjgap may be desirable. This estimation may give 
a guiding principle for materials synthesis. Sequential 
pluses of the longitudinal magnetic fields may remedy a 
quick decay of the sliding. 

The sliding motion may be signaled by the spin-density 
accumulation inside each soliton (kink) and emergence of 
periodic arrays of the induced magnetic dipoles carrying 
the transport spin current.^ From theoretical viewpoints, 
the coherent sliding is always possible to occur as a direct 
consequence of the phase rigidity and Galilean symme- 
try in any type of density waves, including spin/charge 
density waves and even inhomogeneous superconducting 
states. However, in many types of such systems, the slid- 
ing motion does not transport experimentally measured 
quantity.^^ In this respect, it is remarkable that the co- 
herent sliding of the CSL accompanies the dynamically 
generated magnetization. 

Another observable consequence of the sliding is ap- 
pearance of the SMF along the helical axis. We showed 
that the time-dependent sliding velocity Z{t) causes 
time- varying Berry cap which causes the SMF. We 
stressed that the dissipative dynamics plays an essen- 
tial role to drive SMF. A salient feature of the CSL is 
appearance of the strongly amplified SMF which is di- 
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rectly proportional to the macroscopic number of soliton. 
Consequently, the SMF is expected to reach the order of 
mV. As reported in Refj^, CHM and CSL are quite ro- 
bust against structural dislocation and crystal defects. 
Its high stability and robustness are direct manifestation 
of the macroscopic order of spin magnetic moments in 
CHM and CSL state. We hope that the present proposal 
may lead to spintronics application based on chiral mag- 
netic crystals. 
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Appendix A: Static deformation 

Plugging the expressions, 9{z) — Tr/2 + s9{z) and 
ip{z) = (po{z) + sip{z), into the static counterparts of 
Eqs. (|9a|) and (|9b|) . and retaining the first order correc- 
tions with respect to s, we have 



die = {2k^svl^z- 



(2K'^sn^z 



^2) 



(Ala) 
(Alb) 



where the dimensionless variables z — \fP^z/K^ Pz = 
HJJS, (3x = Hx/JS are introduced. Eq. (|Ala|) is 
the homogeneous Lame equation, while Eq. (jAlbp is a 
non-homogeneous Lame equation. The solution of Lame 
equation is well known to be given in a form^ 



6(z) 



(A2) 



where H and Z are Jacobi's eta and zeta functions, re- 
spectively with the parameter a being determined by 
dn^a = -16£;V7r2. 

A solution for the non-homogeneous equation (jAlbp is 
obtained by using homogeneous solutions <^i,2- In the 
inhomogeneous term, we ignore Qp^i as compared with 
Pz- This treatment is justified because Pz and 9i are of 
the same order, and Qq 1. Using a method of vari- 
ation of parameters for non-homogeneous second-order 
differential equation, we readily construct the solution as 

e{z) = y^^'^W-^ (^^2{z) j dz^i{z) - ^i(z) J dz02{z) 

(A3) 

where W is the Wronskian, 



W = ipi{z)ip^{z) - <f2{z)<fi{z). 



(A4) 



into Ea. (IA4[) . we finally obtain 

I" ^1 (#) 



2K 01 



2K J 



Z{a) 



(A7) 



where 6i{x) {i — 1, 2, 3, 4) denote the elliptic Theta func- 
tions and i?4 = 6*4(0) = ^2k' K/tt. Here k is the com- 
plementary modulus. 

Final task is to perform integrals in (IA3|) is performed 
by using the Fourier transformation of (^1,2 (^) to give 
Eq. ((T9)l . The derivation is similar to the calculation of 
Fourier coefficients for the Jacobi sn functionJ^ We start 
with the Fourier transformation 



By definition, we have 
27rc_r 



n— — oo 



(A8) 



\_{x + a) 



' dx. 



To evaluate the integral, we use the contour taken as a 
parallelogram with the corner points — tt, tt, tt + t:t and 
— 27r + 7rr, where t = iK / K. The singular points inside 
the contour are zi = — tt + 7rr/2 and Z2 — ■kt /2. After 
straightforward computations, we obtain 



ci,2n = for even n 

6 

Cl.2n - 



I J£a\ 
4 [2k) 



e[ sinh [tt [nK' T ia) /2K] 
for odd n. Finally we have 



(A9) 

for odd n. (AlO) 



Ci = - 



+ CXD 

n— — 00 



2K 



+ 00 

C2 = E Clr, 



Z{a) 



Z7rn 
2K 



Z{a) ' 



(All) 
(A12) 



where in the summation only terms with odd n are re- 
tained. 
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Appendix B: AC field driven oscillations of the CSL 



We derive a solution of the LLG equations (|9al9bp 
for the periodic longitudinal magnetic field Pz{t) — 
Pzo sin 51t. The solution is sought in the form 6 = tt/2+9i 
and if ~ ipo + ifi^ where the additions 9i, ipi are of the 
same order of magnitude as the magnetic field /J^. To 
provide an analytical treatment, we consider the limit of 
the small when the approximations (Pq{z) ~ Qqz and 
^0 = 7r/2 are relevant. Moreover, we assume a smallness 
of the Gilbert damping, when the problem becomes it- 
erative. At the first stage, we find solutions for 9i, ipi 
at a = and plug them into Eqs. (|9al9bl) to obatin new 
values valid for non-zero a. 

At a = Eas. (|9al9bp read as 



a7 



9Vi 



— — ~<^o''i + 'TTT 



cos QoZLpi, 



fixCosQazl 



(Bl) 



(B2) 



and can be easily resolved through the substitutions 
(^^°^(z,t) = (Ai + A2/3a: cos Qqz) cos rir. 



6*^°^ (2, r) = (Bl + B2px cos Qqz) sin f]r, 

where Ai^2 and B12 are the unknowns. 
This results straightforwardly in 



(B3) 



(z t) - 



(B4) 

cosQo-z sin fir. (B5) 



A requirement of smallness of the corrections amounts to 
/3,o < f!, 2/3,^,oO§ < f^(2Qg-f^'), and /3,/3,o « 2g^- 
il^. Taking Qo ^ 10~^ and J7 ^ 10~^ in dimensionless 
units (or 1 GHz in physical units flro), we can suppose, 
for example, Pzo/^ ~ 0.1, 13^ ^ lO"" (100 Oe), /3^o 
10-5 (10 Oe). 

By assuming smallness of the Gilbert parameter a, we 
organize the iterative procedure to find solutions of the 
system (|9al9bp with the time derivatives in the right- 
hand sides estimated from Eqs. (IB4IB5p . The calculation 
yields 



cos fir (B6) 



1 - 



vi^ - 2g4 



f!2 - 2Q4 



cos QqZ 




FIG. 12: The time-dependent variation of the polar angle 
for the central site (A^ = 50000): numerical data (dots) and 
analytical result (line) given by Ea. (|B6p . The fields are taken 
as in Fig. 8 in the main text, a = 0.1. 
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FIG. 13: The time-dependent variation of the azimuthal angle 
for the central site [N = 50000): numerical data (dots) and 
analytical result (line) given by Ea. (|B7[) . The parameters are 
the same as in the previous Figure. 



4 I cosc</oz 

^^0 



A direct comparison between the numerical results and 
the analytical predictions is done in Figs. [TH [131 Ob- 
viously, there is a good agreement in the steady-state 
regime, when eigenmodes fade away. 



Appendix C: Absence of SMF for rigid motion 

We here note that the dissipationless rigid motion 
of the CSL never causes SMF. The results can be 
obtained from general considerations. In Villain's 
representation^ the spin component — hS cos 9 and 
the angle (p made by the projection of the spin in the 
(x, y) plane are conjugated canonical variables. Eqs. (j9ap 
and (|9bp are written in the new variables acquire the 
Hamiltonian form 



ipi{z,T) — ipf\z,T) + a^^sm.Q,T (B7) 



dip 

at 



'df 



m 

dip ' 



(CI) 
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whereas the fictitious electric field 
follows 



37)1 is presented as 



1 fdS'dip dS'dip 



dt dz dz dt 



(C2) 



The spin motive force generated along the path of length 
L reduces to the contour integral in the phase space of 
the conjugated variables {<f,S^) 



1 



e{t) - / dz Eiz, t)^ — <p ( —dip - -^dS 



d_ 
dip 



2S 
d(p 
'dt 



dip 



^ dS^ \ dt 



dpdS"' 



(C3) 



Plugging Eqs. dCip into the formula we obtain e{t) = 0. 
This rigorous result shows dissipationless Hamiltonian 
dynamics of any spin texture never causes finite SMF. 
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